Low energy fixed points of the cr-r and the 0(3) symmetric Anderson models 
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We study the single channel (compactified) models, the g-t model and the 0(3) symmetric Ander- 
son model, which were introduced by Coleman et al, and Coleman and Schofield, as a simplified 
way to understand the low energy behaviour of the isotropic and anisotropic two channel Kondo 
systems. These models display both Fermi liquid and marginal Fermi liquid behaviour and an un- 
derstanding of the nature of their low energy fixed points may give some general insights into the 
low energy behaviour of other strongly correlated systems. We calculate the excitation spectrum 
at the non-Fermi liquid fixed point of the a-r model using conformal field theory, and show that 
the results are in agreement with those obtained in recent numerical renormalization group (NRG) 
calculations. For the 0(3) Anderson model we find further logarithmic corrections in the weak cou- 
pling perturbation expansion to those obtained in earlier calculations, such that the renormalized 
interaction term now becomes marginally stable rather than marginally unstable. We derive a Ward 
identity and a renormalized form of the perturbation theory that encompasses both the weak and 
strong coupling regimes and show that the x/7 ratio is 8/3 (x is the total susceptibility, spin plus 
isospin), independent of the interaction U and in agreement with the NRG calculations. 
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I. INTRODUCTION 

In developing a simplified way to understand the low 
energy behaviour of the isotropic and anisotropic two 
channel Kondo systems Coleman et al. Q and Coleman 
and Schofield introduced two new single channel im- 
purity models. These are the ct-t model and the 0(3) 
symmetric Anderson model. The cr-r model is similar in 
form to the two channel Kondo model in that a local- 
ized spin s = 1/2 is coupled via an exchange interaction 
to two conduction electron terms, but only one of these 
terms is the conduction electron spin cr, the other term is 
the conduction electron isospin r. The Hamiltonian can 
be expressed in the form, 

H = [Jia(O) + J2f(0)] • 

oc 
n—0.(T 



T-z{n) = ^(4^c„^ + cT„^c„^ - 1), 



(3) 



(1) 



where Ji and J2 are the two exchange couplings, and the 
conduction electrons are in the form of a tight-binding 
chain with cj^^, c^^, the creation and annihilation oper- 
ators at site n, and tn is the nearest neighbour hopping 
matrix element. The spin operators, 

CT+(n) = cl^c^^, a-{n) = c|^c„p 
^zW = ^(ci|C„^ - c|^c„J, (2) 

in the basis spanned by the two singly occupied fermion 
states, |0(t),l(i)) and |l(t),0(|)), give a representation 
for the SU(2) algebra for spin a ~ ^. The isospin opera- 
tors, 

T+(n) = hirci^ci,, r-{n) = (-l)"c„,c„^. 



give a representation of the same algebra in the space 
spanned by the zero and double occupation states, |0, 0) 
and |1, 1). The operators for the two conduction electron 
channels of the two channel Kondo model have been re- 
placed by the operators of a single conduction electron 
channel, so the term 'compactified' has also been used to 
describe the model. The model as introduced by Cole- 
man et al. was for the isotropic case Ji = J2. 

The 0(3) symmetric Anderson model is a modified 
form of the symmetric Anderson model (SAM). The sym- 
metric Anderson model can be written in the form. 



n—0,(7 

+ ^(4aCd. + c^coj + U{nd^ - ^){nai - i), (4) 

where c^^ and c^^ are creation and annihilation opera- 
tors for the localized impurity d state with spin a. The 
electrons in the impurity state interact via a Coulomb 
matrix element U, and V is the hybridization matrix el- 
ement of the impurity state with a conduction electron 
chain which is of the same form as that given in (|l]) . This 
model was modified by Coleman and Schofield by the 
addition of an anomalous hybridization term. 



V,{clc, 



■ Cn I C, 



■ ^i^Oi + 4i '^dil- 



ls) 



The motivation for the addition of such a term is that 
if one applies a Schrieffer- Wolff transformation to the 
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model in the local moment regime (large U) then it maps 
into the cr-r model with Ji and J2 given by 



'J\ — 77 7 'H — 



U 



U 



(6) 



The symmetry of the model can be made explicit by 
the introduction of the Majorana fermion representation. 
The Majorana fermion operators for the d electrons are 
defined by 

•^i = :^(^'iT + 4t) ^2 = -^(cdT - 4t)' 
^3 = ^(cdi+4|) do = ^{c^^-c\), (7) 
which satisfy the commutation relations, 

{c^Q,rf/3} = ^a,/3, (8) 

where {} indicates an anticommutator. 

The Majorana operators for the conduction electrons 
are similarly defined by 



V2 



V2 



-((-f)"c„T+cl^), 
'((-1)Xt-4t)' 



V2 

inrr /2 



((~l)"c„,+cl,) 



le 



V2 



((-i)Xi-ci,), 



(9) 



with commutation relations as in (|^). The modified An- 
derson model in this representation takes the form 

3 

a = l 
3 00 

+ j ini^o^in + l)V'a(n) + Udid2d3do, (10) 

Q=0 n=0 

where Vb = V — 2Vs.. The Hamiltonian is invariant under 
the orthogonal transformation, 

3 3 

da = ^ aa,f3df3, ipa{n) = ^ aa,i3ip p{n) , (11) 
/3=1 /3=1 

where ^a,paa' ,0 = Sa,a' and a, a',/3 = 1,2,3, which 
corresponds to 0(3) symmetry as the Majorana states 
are not included. For Vo = we recover the standard 
symmetric Anderson model (T4 = 0) and the model is 
then symmetric under the corresponding 0(4) transfor- 
mation with the Majorana states included. 

These two models do not conserve either the number 
of particles or the spin. However, there is a combination 



of spin plus isospin, which we denote by jtot, which is 
conserved. For the Anderson model this corresponds to 



jtot = Sd + fd + 



(12) 



n=0 



where fd is the isopin operator for the impurity site. As 
the a-T model corresponds to the localized limit in which 
the impurity charge fluctuations are suppressed the fd is 
omitted for this model. Clearly this combination satisfies 
the usual SU(2) angular momentum commutation rela- 
tions, so the many-body states of these models can be 
classified by the quantum numbers jtot and jtot,z- The 
combination of spin and isospin for a particular site jn 
when expressed in terms of the Majorana fermions takes 
the form, 



(13) 



where ijj{n) — (jjji{n)^ip2in),'ip3{n)). Hence we will re- 
fer to the a — 1,2,3 Majorana components as the com- 
ponents of a vector field, and the a = as the scalar 
component. 

The a-T model differs from the two channel Kondo 
model in that there is no possibility of overscreening the 
impurity spin. Overscreening occurs in the two channel 
case in the strong coupling limit Ji — > cxd, J2 cxd, 
when a conduction electron in each of the two channels 
binds to the impurity spin resulting in a localized state at 
the impurity which again has spin 1/2. In the Ji — s- 00, 
J2 00 limit of the cr-r model, however, because the spin 
and isospin configurations are mutually exclusive states 
of the same channel the impurity cannot bind to both 
simultaneously and overcompensation cannot occur. 

Coleman et al. established a correspondence between 
the two models in perturbation theory (after bosoniza- 
tion) and they conjecture that the low energy physics of 
the two models is essentially the same (see also ||] ) , with 
the extra states of the two channel model playing no sig- 
nificant role. For the isotropic case, Ji = J2 = J they 
obtained an effective Hamiltonian for the cr-r model by 
assuming that the low energy fixed point corresponds to 
J = 00, and then calculating the leading order terms 
in an expansion in D/J where D is the band width 
of the conduction electrons. The effective Hamiltonian 
to leading order in D/J has an interaction term which 
corresponds to that in the 0(3) Anderson model with 
U ~ / and Vb = 0, so that corrections to the fixed 
point correspond to a low order perturbation expansion 
in powers of U for the 0(3) Anderson model. The sec- 
ond order terms generate InT terms in the susceptibility 
and specific heat, and a Wilson ratio for the logarithmic 
terms of 8/3, which is the low energy behaviour of the 
two channel Kondo model. 

There are several reasons why we think these models 
are worthy of further study. The low order perturbation 
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theory in U (for Vq = 0) gives a self-energy corresponding 
to the marginal Fermi liquid form whose imaginary part 
is given by 



ImI](2)(cj,T) = - 




and real part, 

ReE(^Hc.,r)«(-^)^ln(J), (15) 

where x = max(|ci;|,r), A — TrV'^pdO) and Pc{0) is the 
density of conduction electron states at the Fermi level 
ii. 

Marginal Fermi liquid theory was put forward in a phe- 
nomenological interpretation of the behaviour of the nor- 
mal state of the high Tc superconductors so it is of some 
interest to understand fully a microscopic model which 
displays this behaviour |^. In higher order perturba- 
tion theory diagrams were found that give logarithmic 
contributions to the irreducible four vertex (at zero fre- 
quency) with a sign indicating that the weak coupling 
fixed point 0| might be unstable. It is important to 
understand the role of these terms in the low energy 
behaviour of the model, and to determine whether or 
not the marginal Fermi liquid behaviour corresponds to 
a stable fixed point. Apart from the non- Fermi liquid be- 
haviour of the isotropic model ( Ji — J2) it has been also 
conjectured, based on Bethe ansatz calculations, that the 
behaviour of the anisotropic two channel model also dif- 
fers from that for a Fermi liquid fixed point Q . If the a-r 
model and the two channel Kondo model have the same 
low energy fixed point then the anisotropic a-r should 
behave in a similar way. These are points that need fur- 
ther clarification. 

In the derivation of the low energy behaviour of the a-r 
model based on the 1/J expansion there are points that 
are not clear. For the single channel Kondo model it has 
been shown by Nozieres [|7|] that, though the low energy 
fixed point corresponds to J = cxi, the low energy physics 
cannot be calculated from a 1/ J expansion (which gives 
a non universal Wilson ratio). We know that for this 
model there is only one energy scale which is the Kondo 
temperature Tk, and so the Wilson ratio (i?) has a uni- 
versal value (i? = 2). The effective Hamiltonian that 
determines the low energy impurity behaviour depends 
only on Tk, and asymptotic expansions about the fixed 
point are in powers of T/Tk (or frequency or magnetic 
field compared to the Kondo temperature). The leading 
interaction terms in the effective Hamiltonian generated 
by a 1/J expansion do satisfy the Wilson criteria for the 
leading corrections to the fixed point because (a) they 
involve operators local to the impurity, and (b) they are 
consistent with the symmetries of the model. However 
the combination of these terms generated in this effec- 
tive Hamiltonian (D/J <C 1), does not correspond to the 



correct low energy effective Hamiltonian for the Kondo 
model. The 1/J effective Hamiltonian, as to be expected, 
reproduces the £)/ J <C 1 physics of the model with a non- 
universal Wilson ratio dependent upon D/ J, and so does 
not give the correct low energy behaviour about the fixed 
point of the weak coupling model. The applicability of 
the 1/J expansion to derive the excitations about the 
fixed point of the weak coupling ct-t model needs to be 
clarified. 

In this paper we bring together information on the low 
energy behaviour of these models, obtained from differ- 
ent many-body techniques, conformal field theory, nu- 
merical renormalization group and perturbation theory 
in order to clarify and resolve some of these outstanding 
issues. We first of all consider the nature of the exci- 
tations in a free single component Majorana chain with 
various types of boundary conditions, and then calcu- 
late possible forms of the many-body excitation spectra 
when the component systems are put together. We re- 
late these results to those obtained using conformal field 
theory. The conformal field approach can also be used 
to calculate the many-body excitations of the isotropic 
a-T model for a particular value of the exchange inter- 
action J — J* ^ corresponding to strong couphng. We 
show that the excitations of the model at this value of 
J correspond to those found at the fixed point in recent 
numerical renormalization group calculations ||^. From 
the conformal field theory we also obtain the form of the 
operators which give the leading order corrections to this 
fixed point. 

For the 0(3) Anderson model we find additional log- 
arithmic terms to those found in earlier work [Q, and 
when these are included the renormalized interaction is 
marginally irrelevant rather than marginally relevant so 
the marginal fixed point is stable. We also derive a Ward 
identity from which we can deduce the Wilson ratio for 
the total susceptibility (spin plus isospin) of the impu- 
rity to the specific heat coefficient. We show that a 
renormalized perturbation expansion for the Fermi liquid 
fixed point (Vo ^ 0) can be generalized for the marginal 
Fermi liquid fixed point (Vq = 0). The renormalized per- 
turbation theory encompasses both the strong and weak 
coupling limits. Finally we relate our results to those ob- 
tained by other methods, and to the results for the two 
channel Kondo model. 

II. FREE MAJORANA FERMION CHAIN 



A Fermi liquid fixed point of an interacting system is 
characterized by quasiparticle excitations which are in 
one-to-one correspondence with the single particle ex- 
citations of the non-interacting system. In the models 
we are considering here the single particle excitations 
of the non-interacting system are independent Majorana 
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fermion excitations. As a preliminary to investigating 
the fixed point of the interacting systems we look first of 
all at the possible excitations of a free Majorana chain. 
We consider just a one component (a = 0) chain with N 
sites described by the Hamiltonian, 



Af-l 



(16) 



n=0 



where we have taken the hopping matrix-element t (real) 
to be the same between all neighbouring sites. The num- 
ber of sites is assumed to be even and we use periodic or 
antiperiodic boundary conditions ipoi^N) = ±V'o(0). 
The Fourier-transformed operators are defined as 



N-l 



e^«'="Vo(n), 



(17) 



We consider first the case of periodic boundary condi- 
tions. In this case k takes only the integer values, 



N 



N 

-1,0,1,...,-. 



(18) 



The Hamiltonian (|T^) can be diagonalized by using Eq. 
H^Yl ite'^'MkM-k) + ^ ite''^''Mk)M-k) 

k<0 k>Q 



- N - N 



(19) 



Due to '00 ("y) — ■'/'o("f"), the last two terms cancel. 
There is a Majorana fermion state at the Fermi level but 
it does not appear explicitly in the Hamiltonian as we 
take 11 = 0. For fc > we define the operators 



(20) 



Using Eq. (|l^ and the commutation relation for the Ma- 
jorana fermions ■^/'o(?^), one obtains 



(21) 



fc>0 



(24) 



In the limit of large N, the single particle excitation en- 
ergies Ck are given by 



£fc — ^ J ^' ^ ~ 1:2,3,, 



(25) 



for small k/N. There is a similar linear branch for k/N ^ 
^ . The former we will refer to as the left branch and the 
latter the right branch. Here we have expressed t in terms 
of the Fermi velocity v-p = 2t and replaced the number 
of sites N by the length of the chain L = N = 21 (we 
introduce I to be consistent with the notation used in 
Sec. HI). 

The ground state energy is given by 



sin(f) 



v-pl 



(26) 



In the case of antiperiodic boundary conditions, k can 
only take the half-integer values 



N 1 
T+2' 



1 1 

2' 2' 



N 



(27) 



In contrast to Eq. (|19|) there is no degree of freedom at 
the Fermi level, so the Hamiltonian reads 



H 



{N-l)/2 

X! <^ki^l{k^)i^f){k) 

k=l/2 



E, 



(28) 



with the dispersion Eq. ( p3| ) and the ground state energy 
-ztVe'^fc. (29) 



k>0 



In the limit of large N, the single particle excitation en- 
ergies are given by 



ttwf , , 13 5 
^k^—k, fc = 

For the ground state energy we find 



(30) 



which is just the anticommutation relation for ordinary 
fermion operators. That means that we can treat the 
operators "00 (fc) and V'o(fc) for fc > as fermion operators 
in the following discussion. The Hamiltonian now takes 
the simple form 



H = J2^ki'lik)Mk)+Eg^p, 



k>0 

with the dispersion 

„ . 27r, 
€k = 2t sm — k 

N 

and the ground state energy 



(22) 



(23) 



E =-t ^ 

sm(f) TT 



1 t:' 
'+2AP 



(31) 



The difference between the ground state energies for pe- 
riodic and antiperiodic boundary conditions is 



E„ 



Er, 



1 

12 



1 

24 



1 nvp 
8~' 



(32) 



where both linearized branches of ( |2^ ) have been taken 
into account. With the restriction to only one branch of 
Majorana fermions, the single particle spectrum is given 
by Fig. 1(a) and 1(b) for periodic and antiperiodic bound- 
ary conditions. 
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From the single particle spectrum we can also build up 
the spectrum for the many-body excitations (relative to 
the ground state) of the scalar Majorana fermion chain, 
which corresponds to progressively creating 'particles' in 
the states with k> 0. The ground state is a vacuum state 
with no particles. The 'particles' which are created are 
combinations of the usual particle and hole excitations of 
a free electron system, which are degenerate as a result 
of the particle-hole symmetry. The resulting many-body 
energies -Eox and their degeneracies dg are given in Tables 
I and II, where Uk denotes the occupation of the single 
particle state of energy e^. 

We can build up the many-body excitations of the vec- 
tor Majorana fermion chain corresponding to the Hamil- 
tonian. 



N-l 3 



= it ^ ^ ipa{n)ipa{n + 1), 



(33) 



in a similar way. The single particle spectra for this 
model for periodic and antiperiodic boundary conditions 
are given in Fig. 2(a) and 2(b), respectively. The number 
of excited states per energy level is simply multiplied by 
three as compared to Fig. 1. The resulting many-body 
states and their degeneracies are given in Table III and 
IV. 



III. CONFORMAL FIELD THEORY (CFT) 

In this section, we apply a current algebra approach 
to the channel isotropic <t-t model along the lines devel- 
oped by Affleck and Ludwig |^ for the two channel Kondo 
model. In our case, however, a new separation scheme is 
used in which the conduction electrons, expressed as Ma- 
jorana fermions, are separated into a density correspond- 
ing to total spin (spin plus isospin, the vector channel) 
and a single (scalar) component Majorana fermion. As 
in the other examples where this technique is used the 
model will be solved at a particular finite but large value 
of the coupling J*. At this value of the coupling the 
impurity spin can be absorbed into the conduction elec- 
trons, and the algebra which determines the excitation 
spectrum is then of the same form as that for the non- 
interacting system. 

The (T-r Hamiltonian can be written as 



H = it^^ tpa(n + l)^a(n) + J jo 



Sd- 



where j„ is the combination of spin and isospin as defined 
earlier in Eq. (|l^). In this form the vector and scalar 
Majorana fermion terms are decoupled so that the model 
can be divided into two parts (see Fig. 3) 



Hsc = it'^ipo{n + l)V'o(n), 

n=0 
oo 3 

^^vcc = it^^ipa{n+ l)i>ain) + J jo ■ Sd- (34) 



If the lattice spacing is a we can define a spatial coor- 
dinate X = na for each spin and in the limit as a — > 
we use Eqs. (|^) and (|^) to define both continuum spin 
and isospin densities with components, a"'{x) and t'^{x), 
where 7 = x,y,z. These densities satisfy a level k = 1 
SU{2) Kac-Moody algebra: 



s%x),s^{x')] 



ik 



= ie°''^-^ {x)6{x - x') + —5c.i3S'{x - x'). (35) 

where s^{x) = cr^(a;), r^(a;), and 5' is the derivative of the 
delta function. If we define similarly a density j{x) for 
the combined spin and isospin then the three components 
of this density satisfy an SU{2) level k = 2 Kac-Moody 
algebra. 

In the continuum limit, the linearization of the dis- 
persion relation can be made, and the conduction elec- 
trons can be expressed in terms of left- and right-moving 
Majorana fermions satisfying the conditions ipa,L{Q) = 
'0Q.i?.(O) with a = 0, 1, 2, 3. In fact, one can continue the 
Majorana fermions to the negative cc-direction by defin- 
ing il)a,L(x) — iJJa.Ri^x), SO the Hamiltonians can now be 
written in terms of left-moving fermions only on the full 
X-axis. The Hamiltonian that describes the excitations 
for the scalar part becomes 

/oo 
dx : i}}o{x)idxtpQ{x) :, (36) 
^00 

where the subscripts "L" have been omitted, the Fermi 
velocity is v-p = 2ta, and : : indicates normal ordering. 
This Hamiltonian has the same form as the one which 
describes the two dimensional Ising model at the critical 
point 0. 

For the vector part it is possible to apply the stan- 
dard technique of conformal field theory and write the 
non-interacting part in the Sugawara form |jll| , which is 
quadratic in the densities rather than the fields, so that 
-ffvec takes the form, 



ZTT 



dx 



i : jix) ■ j{x) 



''^^6{x)j{x) ■ Sd 



(37) 



which is normal ordered with respect to the non-inter- 
acting system. This Hamiltonian is expressed entirely in 
terms of the combined spin density j{x) and the impurity 
spin, and is the only part which contains the interaction 
term. 
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We choose boundary conditions V'a(0 = e V'q(— 0' 
on a large circle —l<x<l for a = 0,1,2,3. The 
momenta are then given by km = j {m+ |:) and the 
energy levels with respect to the chemical potential arc 
expressed by Em = v^km with m = 0,±1,±2,- • •. In 
the case of particle-hole symmetry, 6 = corresponds to 
the case of odd number of sites (m € Z, integers) in the 
lattice version, while 6 = tt/2 to the case of even number 
sites (to G Z + i, half-odd integers), these correspond 
to periodic boundary condition (PBC) and anti-periodic 
boundary condition (APBC), respectively. 

We now look at the J = case with APBC and PBC. 
After taking a Fourier transform, the scalar part becomes 



H. 



(0) 



E 

m>0 

TTVp Cue 

I 24 ' 
I 12 ' 



TO e Z 
TO e Z 



(38) 



where Csc = 5 is the central charge of the Virasoro alge- 
bra, 

[Ln, Lra] = [u - to)L„+™ -|- ^^("-^ " n)6n+m,0, (39) 

where Ln is the nth component of the mode expansion of 
the stress-energy tensor, which can be used to generate 
the excitations of the system (for details see [|2|). The 
eigenvalues of this Hamiltonian can be expressed in the 
form 



Esc — 



48 ^ ^"'"s 



where Ajg 



0, 



ising — ") 5, jq- This notation is that used in 
the context of the two dimensional Ising model where 
the values of Aisi„g correspond to the scaling dimensions 
of the primary fields: identity operator /, energy density 
operator e, and Ising order parameter a, respectively. For 
APBC, the first two primary fields appear and each with 
degeneracy dgsc = 1. In the PBC case, only a appears 
with a degeneracy dgsc — 2. These results agree with 
those obtained in Tables I and II in the previous section. 
Those in Table I correspond to Ajging = -j^, and in Tabic 
II correspond to Ajsing = 0, i. 

In terms of the Fourier transform of the spin plus 
isospin density. 



where Cvec = § is the central charge for the vector field. 
The eigenvalues for this part of the Hamiltonian which 
describes the vector Majorana fermions can be written 
in the form, 



Ev. 



TTVp 



1 

-16 + 



J-(J + 1) 



where j is a good quantum number and tt-voc is a non- 
negative integer. The primary fields of the conformal 
field theory corresponding to this model are j = (sin- 
glet), j = 5 (doublet), and j = 1 (triplet), with scahng 
dimension j{i + l)/4 and degeneracy dgvoc = (2j + 1). 
In the PBC case j — \ and the results correspond to 
those given in Table III for the excitations relative to the 
ground state, while for APBC, j = 0, 1 and they corre- 
spond to those given in Table IV. The 0(i) corrections 
to the ground state energy, which depend on the central 
charges Csc and Cvec, are also in agreement with those 
found in the previous section. 

When we put the vector and scalar parts together, the 
finite-size excitation spectrum of the full model is given 
by 



Aising) — Esc + Eycc 
J(J + 1) 



TTV-p 



1 

12 



^ Ising 



(41) 



Not all the combinations of the spin and charge degrees 
of freedom are allowed. Only those combinations rep- 
resenting composite fermions corresponding to the 0(4) 
symmetry of our original free fermion model are allowed. 
To satisfy this condition the vector part and scalar part 
should have the same boundary conditions. Therefore, 
the possible combinations of (j, Ajsing) are (0,/), (0, e) 
or (1,/), and (l,e) for the APBC (to G Z + \), while 
(i,(T) for the PBC (to G Z). The finite-size excitation 
spectra of the non-interacting case for PBC and APBC 
are shown in Table V and Table VI separately. 

We now turn to the case of the interacting model J ^ 
0. For a particular value of the coupling J = J* = vy 
it is possible to incorporate the impurity spin and define 
a new spin density j'{x) = j{x) -\- 2'KS(x)sd such that the 
impurity spin operator Sd disappears explicitly from the 
Hamiltonian. The vector part of the Hamiltonian is now 
formally the same as for the non-interacting model. 



1 

2^ 



m— — oo 



the free vector part may be cast into the form, 



^vcc 2^ ^ ^ j ( 

m—l 



TTVf 1 



■ j{m) 



-^%f, TOGZ + l 
TOGZ 



(40) 



The new spin-density operator j'{x) obeys the same 
SU{2)2 commutation relations as the old one i{x) so that 
at the special point J = J* the same spectrum of excita- 
tions are generated for Hsc + H*^^ in the conformal field 
algebra as in the non-interacting case. The value of the 
coupling J* is very large, being of the order of the band 
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width of the conduction electrons, so it corresponds to a 
strong couphng hmit of the original model (see Fig. 4). 
We are interested in the weak coupling limit of the orig- 
inal model. However, it might describe the fixed point 
Hamiltonian in a Wilson type of renormalization group 
calculation (where the energy scale is progressively re- 
duced by eliminating the higher energy states), similar 
to the case of the s-d (Kondo) model where the fixed 
point is the strong coupling one. This was conjectured in 
earlier work ^ , and would be consistent with the expec- 
tation that the fixed point Hamiltonian should be con- 
formally invariant. In the next section we look at the 
results of explicit numerical renormalization group cal- 
culations that confirm this. Before doing so, however, we 
will look at physically relevant excitations of the model 
at J = J* . These differ from those already considered for 
the non-interacting system because, when the impurity 
spin has been absorbed, the allowed combinations of the 
vector and scalar excitations are not the same. When 
we change from j{x) to j'{x), an extra spin ^ degree of 
freedom is added to the SU{2)2 spin density whereas the 
quantum numbers of the scalar part are unchanged. The 
changes in the vector part follow the SU {2)2 fusion rule 

As a result, the boundary conditions of the vector part 
are changed compared to the non-interacting cases: the 
APBC changes to PBC, while PBC to APBC. Therefore, 
the boundary conditions of the vector chain and scalar 
chain are no longer the same. We now have two sets 
of excitation spectra, which correspond to the combina- 
tion of the vector part with PBC and scalar part with 
APBC, or vice versa. These two spectra, which are dis- 
tinguishable, are displayed in the Table VH and Table 
VHI, respectively. We will relate these results to those of 
the numerical renormalization group in the next section. 

IV. NUMERICAL RENORMALIZATION GROUP 
(NRG) 

The numerical renormalization group method was de- 
veloped by Wilson for the Kondo problem and later 
applied to the single impurity Anderson model by Kr- 
ishnamurthy et al. [Q. It is based on the mapping of 
the model to a semi-infinite chain with the spin or im- 
purity on one end of the chain coupled to the first con- 
duction electron site via an exchange coupling J or a hy- 
bridization V . The hopping matrix elements i„ between 
neighbouring sites along the chain approach a constant 
value for large n as long as the mapping is performed ex- 
actly. In the numerical renormalization group, however, 
certain irrelevant degrees of freedom are neglected which 
results in an exponential decay f„ (x A~"/^. The origin 



of the parameter A is the logarithmic discretization of 
the continuous conduction band with the energy mesh 
e„ = ±A^", n = 0, 1, 2, . . .. Details of the derivation of 
the semi-infinite chain form used here can be found in 
iil and 0. 

Important for the discussion in this paper is the A- 
dependence of the results due to the logarithmic dis- 
cretization. As has been shown earlier, the physical prop- 
erties calculated from the NRG depend only weakly on 
the discretization parameter A as long as the number of 
states taken into account is not too small (usually 500 - 
1000). On the other hand, the many-body energy levels 
at the fixed points are A-dependent. In order to com- 
pare the fixed-point spectra with the finite size spectra 
obtained from the CFT, we have to take the limit A 1. 
This will be discussed in the following subsection. 

A. The non-Fermi liquid fixed point 

We want to concentrate on the fixed point spectra it- 
self, not on the deviations responsible for the InT-terms 
in the specific heat and the susceptibility which will be 
considered in Sec. V. A flow diagram (see also ^) show- 
ing the lowest lying many-body energy levels is presented 
in Fig. 5 for the isotropic case T4 = V/2. 

The value of the discretization parameter is A = 2.5. 
This diagram shows the crossover from the free-orbital 
fixed point via the local-moment one to the non-Fermi 
liquid fixed point. The influence of the local moment 
fixed point is not great due to the relatively small value of 
U . The excitations of this fixed point are quite different 
to those found at the Fermi liquid fixed point for 14 7^ 
V/2. 

In the previous section, the ct-t model was written 
in terms of Majorana fermions and it was pointed out 
that in the isotropic case the scalar part of the conduc- 
tion electron chain is decoupled from the impurity, and 
that the strong coupling affects only the vector part. As 
J ~ V^^/J7, this corresponds to the limit V ^ 00 oi the 
Anderson model, where V is the hybridization between 
the vector part of the impurity and the first conduction 
electron site (Vb = 0). If we start with a total odd num- 
ber of sites (including the impurity) the effective number 
of sites remaining is odd in the vector part and even in 
the scalar part (see also Fig. 4). The effective reduction 
in the length of the conduction chain at strong coupling 
changes the excitation spectrum of the vector part in the 
same way as a change of boundary condition on the free 
periodic chain. Hence, the combined many-body excita- 
tions at this strong coupling point are the same as those 
calculated via the CFT (see Table VIII). 

As the model used in the NRG calculations depends 
on A through the hopping matrix element i„ ^ A^"/^, 
the single particle excitations for the free chain for this 
model depend on A. Our calculations for the free chain 
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in Sec. II and III, however, correspond to A = 1. Having 
identified the low energy fixed point of the a-r model in 
terms of the free Majorana chain (with changed bound- 
ary conditions), we can use this information to build up 
the many-body excitations for A — 2.5, from the single 
particle ones (A = 2.5) and compare these with those 
found in the numerical renormalization group calcula- 
tions {N — > oo). In Fig. 6, the lowest lying many-body 
energies are calculated from the single particle ones for 
A in the range, 1 < A < 4. We can see from this figure 
that the non-equally spaced excitations go over smoothly 
to the equally spaced spectrum as A 1. The complete 
agreement of this spectrum (calculated from the single 
particle spectrum A = 2.5) with that calculated from the 
full NRG calculation confirms the identification of the 
low energy fixed point of the cr-r model with the strong 
coupling limit. (However, the NRG gives a ground state 
degeneracy of four in contrast to the CFT result dg — 2 
(see Table VIII); the origin of this discrepancy is still 
unclear.) 

We want to point out that, in our case, in order to 
obtain Fig. 6 it was not necessary to run the full NRG- 
program in the entire A-range. As soon as the relation of 
single particle and many-body energies has been verified 
for a certain value of A the many-body energies can be 
calculated from the single particle ones for A = 1 to 
compare with the CFT. This is in contrast to the work 
of Affleck et al. who studied the two channel Kondo 
model with the NRG. In their case, no relation was made 
between the fixed point spectrum and the single particle 
energies of free conduction electron chains so to find the 
many-body energies to compare with the CFT they had 
to run the full NRG program for A — > 1. This causes 
numerical problems as the sequence of iterations does 
not converge in this limit. 

If we start off with a total even number of sites, we get 
agreement with the CFT results with periodic boundary 
conditions for the vector part and antiperiodic boundary 
conditions for the scalar part (see Table VII). 

B. The Fermi liquid fixed point 

For any value of T4 7^ V/2 the system flows to the 
Fermi liquid fixed point of the standard single impurity 
Anderson model defined by T4 = (see the discussion in 

) . Both scalar and vector part of the first conduction 
electron site are now strongly coupled to the impurity so 
that the number of effective sites is the same in both the 
remaining scalar and vector part. The combination of 
scalar and vector part with the same boundary condition 
describes the Fermi liquid fixed point, as discussed in the 
previous section. Again, both the fixed point spectrum 
and the degeneracies (in the limit A ^ 1) are in agree- 
ment with the CFT results given in Tables V and VI. 



V. LEADING ORDER CORRECTIONS TO THE 
FIXED POINTS 

We are primarily interested in the behaviour of the 
models at the non- Fermi liquid fixed point, which corre- 
sponds to J I — Ji for the cr-r model and 14 = y/2 for 
the 0(3) Anderson model. Low order perturbation the- 
ory and the numerical renormalization group calculations 
give a zero point entropy of | ln2, a logarithmic temper- 
ature dependence of the specific heat C(T) oc TlnT , 
and a logarithmic temperature dependence of the 'spin 
plus isospin' susceptibility Xt^) InT. In this section 
we look at arguments that help us to identify the leading 
irrelevant operators at this fixed point that generate this 
behaviour. 

Because the low energy fixed point of the cr-r model 
is the strong coupling one, as is the case for the one 
channel s-d or Kondo model, we can use the analo- 
gies between these two to derive an expression for the 
leading irrelevant operator responsible for the low tem- 
perature behaviour. For the ordinary single-channel 
Kondo model, the leading irrelevant operator is sim- 
ply J^{x — 0), the squared conserved spin density op- 
erator, which is the only SU(2) spin invariant opera- 
tor with scaling dimension 2. For the present reduced 
model, we know that the leading irrelevant operator al- 
ready exists in the non- interacting case where i?voc, in 
terms of three-component Majorana fermions with an 
£'[7(2)2 Kac-Moody algebra, has a supersymmetry as- 
sociated with the invariance under an exchange of the 
spin and isospin degrees of freedom, and the correspond- 
ing supersymmetric current operator can be defined by 
Js{x) = i^/2ip{x) + Oj{x), where 6* is a coordinate which 
anticommutes with the Majorana operators (0^ = 1) 
[ p^ . Thus, we have (x) 9'4'i{x)4>2{x)4>3{x), and from 
the free field theory we know J^{x) is the unique 0(3) 
invariant operator with scaling dimension 3/2. There- 
fore, we are led to identify the leading irrelevant oper- 
ator for the reduced single-channel impurity model as 
J^(a; = 0) ~ 6'V'i(0)i/'2(0)V'3(0). As is an anticom- 
muting coordinate one can make a connection with the 
leading irrelevant operator as calculated by Coleman et 
al. §. In the 1/J expansion they find a leading order 
interaction term of the form $(0)'(/'i(l)V'2(l)''/'3(l) with 
$(0) = ■!/'i(0)V'2(0)?/'3(0). As $(0) is an anticommuting 
operator it acts in the same way as 9, so one can conclude 
that they are essentially the same. We discuss this more 
fully later. 

The numerical renormalization group results we have 
described are for the 0(3) Anderson model which, as 
we discussed in the introduction, is equivalent to the 
CT-r model in the large U (local moment) regime with 
J ~ V'^ /U. We have, therefore, interpreted the strong 
coupling fixed point, in Anderson model terms, as cor- 
responding to the limit V — > 00, which effectively re- 
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moves the first conduction site from the vector Majorana 
fermion chain. We end up for both models with the fixed 
point described by free vector and scalar chains of Majo- 
rana fermions with effective lengths which differ by one 
site. However, the NRG results for the Anderson model 
for V^ = V jl show the same non- Fermi liquid behaviour 
/or all values of U. It is possible to give an equivalent 
interpretation of this fixed point as the U = limit of the 
Anderson model. This situation is essentially the same 
as the V ^ oo one as it corresponds to free vector and 
scalar Majorana fermion chains whose lengths differ by 
one site, as only the vector Majorana chain is coupled to 
the impurity. A similar situation occurs for the standard 
single channel Anderson model where the fixed point can 
be interpretated as a Fermi liquid fixed point for all val- 
ues of U. We develop this approach further in Sec. VII. 

We now briefly look at the arguments used in the NRG 
to identify the leading irrelevant operator at the non- 
Fermi liquid fixed point for the 0(3) symmetric Anderson 
model. We take Hn to be the non-interacting model 
{U = 0) with an N site conduction chain, and consider a 
perturbation of the form, 

dHi = -;7^o(-l)^i(-l)^2(-l)V^3(-l) (42) 

where V'o(~l) is an additional Majorana-fermion that 
does not couple to the scalar chain starting at site (see 
Fig. 7), so our effective impurity site here corresponds to 
= — 1. Wilson devised a way of determining the de- 
gree of irrelevancy (or relevancy) of a particular operator 
by considering how the operator scales on increasing the 
length of the chain N N + 2, when expressed in terms 
of the creation and annihilation operators for diagonal- 
ized states of H^^ and respectively. Because of the 

overall scaling factor the appropriate terms to consider 
are SH' = A^dHi{N) and SH" = A^6Hi{N + 2). 
Following this line of reasoning we find 

SH^ ^ A-^^^dH[, (43) 

so that the effective perturbation on the fixed-point 
Hamiltonian is reduced by a factor A~^/^ when N is in- 
creased by 2. This result is in contrast to the standard 
SAM where a perturbation of the form ( ^2| ) decreases by 
a factor of A~^. This difference is due to the coupling of 
all four components of the Majorana fermions to the rest 
of the chain. Another difference is that (^) is the only 
leading irrelevant perturbation to the non-Fermi liquid 
fixed point whereas in the standard case a perturbation 
in form of a hybridization between site -1 and site also 
reduces with A~^. 

In the 0(3) symmetric Anderson model the hybridiza- 
tion term that couples the vector parts of site -1 and 
has the form 

3 

5H2=iVY,M-'^)M0)- (44) 



Using the same type of argument as just given we can 
show that the corresponding scaled quantities, 5H2, SH!^, 
for the N and -I- 2 site models are related via 

5H'^ = K-^5H'^ . (45) 

The perturbation 5H2 therefore is more irrelevant than 
5Hi . We know from general arguments advanced by Wil- 
son that other local operators involving more sites, or 
more complicated operators at the initial sites, only give 
terms with a higher degree of irrelevancy. 

In the definition of 5H2 we did not take into account 
the coupling between the scalar parts of site -1 and 

5H-i = iVaM-^)M^) ■ (46) 

This perturbation scales as SHf^ — A^^^SH^ and drives 
the system away from the non- Fermi liquid fixed point. 
However, it only exists in the anisotropic case 14 7^ V/2 
where we know from the NRG results that the non- Fermi 
liquid fixed point is unstable. 

Here we have regarded the fixed point to correspond to 
the U = model. If we had considered it to correspond 
to y ^ 00, nothing would have essentially changed in 
the argument except the way we label the initial sites in 
the chains. 

Putting together these leading irrelevant terms with 
the fixed point Hamiltonian gives us a renormalized form 
of the 0(3) Anderson model. We postpone discussion of 
how to parameterize this model, and how to calculate the 
low temperature thermodynamics about the fixed point 
to Sec. VII. Before that we consider some results for the 
0(3) model is its "bare" form as defined in Eq. (|To|). 

VI. PERTURBATION THEORY AND WARD 
IDENTITIES 

In this section we consider the perturbation theory in 
powers of U for the 0(3) Anderson model. In an ear- 
lier paper it was shown that there are diagrams for 
the model with Vq = that give vertex corrections that 
diverge as ln(T) as T ^ 0, indicating that the U term be- 
haves like a marginally relevant operator. Here we show 
that there are other vertex corrections of the same order 
which also give ln(T) terms, such that the net effect is 
that the ln(r) terms cancel. This alters our conclusions 
about the stability of the fixed point in the presence of U. 
We show that similar terms that lead to higher powers of 
ln(r) for the susceptibility Xj also cancel. We also derive 
a Ward identity that allows us to demonstrate that the 
Wilson ratio for Xj/l isf independent of U, in agreement 
with the NRG resufis ||. 

The most direct way to exploit the symmmetries of the 
model, and deduce a Ward identity, is to use the func- 
tional integral approach, which can also be used to gener- 
ate the perturbation theory. We can make the standard 
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transformations and express the partition function Z for 
the model as a functional integral over Grassmann vari- 
ables associated with both the impurity and conduction 
electron Majorana states, 



/3 CXD 
^ n „ n 



(47) 



a=0 n=0 

where the action S is given by 



In the 0(3) model we have Vi = V2 = V3 = F so 
we denote Aq = A for a = 1,2,3. When Vq = in 
Eq. (|l^) the localized d, a = Majorana fermion is de- 
coupled from the conduction electrons and Aq — 0. In 
the functional integral (^) the a = Majorana fermion 
then has a zero energy mode at T = and scattering of 
the other Majorana fermions with this zero energy ex- 
citation leads to infrared singularities and to non-Fermi 
liquid behaviour ||l|,H. 



S 



r 1 

CXD 

+ ^ i^ain, T)dTipa{n, r)} 

n=0 

3 00 

+ ^ ^ tn1pa{n + 1, T)lpa(n, t) 



a=0 n=0 
3 

+iV ^a(0, r)d^{r) + «K)V'o(0, T)do(T) 



A. Low order perturbation theory 

We can develop the perturbation theory in U of refer- 
ence [H for the impurity Green's functions in the func- 
tional integral formalism from the generating function 
Z{rja(T)). This is constructed by adding a coupling to 
a local Majorana Grassmann field r]a{T) to the reduced 
action (p9|) of the form. 



Q = l 



+ C/dl(T)d2(T)d3(T)do(T) 



(48) 



/ dr rja {T)da{T). 
•^0 a=0 



(51) 



Because the action is purely bilinear in the conduction 
electron Grassmann variables these can be integrated 
over to give a reduced action S^ed which is expressed in 
terms of the impurity Grassmann variables only, 



'S'rod — 



I f dT f dT'y^d^{T)[dr5{T^T')+RdT-T')]d^{T') 
^ Jo Jq n 

(49) 



The Green's functions are then obtained by taking the 
appropriate functional derivatives with respect to the 
fields. The thermal Green's function G^°j(r — r') for the 
impurity Majorana fermions for [/ = is given by 

GS(T-r') = -(T.d.(r)dM^')>o 

= 4 y^Gi''j(c.„)e— (52) 



where 



+ U / drdl{T)d2iT)d3{T)do{T). 

Jo 

This expression differs from that for an isolated impu- 
rity only through the term Ra{T — t'). The form for 
Ra{T — t') can most conveniently be expressed in terms 
of its Fourier coefficient Ra{uJn)- The functional integral 
can be re-expressed as an integration over the Fourier 
coefficients of the Grassmann variables dct(w„), where 
dcir) = l//3^„d,(c^„)e-*""" with LUn = (2n + l)7r//3 
so as to satisfy antiperiodic boundary conditions d^ (/3) ~ 
—da{0). The Fourier series for Rair — r') is Ra{T — T') = 
l/Z^Z^n ^ct('^«)'^~"^"'^~^ ^. For a conduction band po{e) 
with —D < e < D, Ra{uj„) is given by 



GS(-n) 



iujn + iAQSgn(iwn) ' 



(53) 



with ujn = {2n + l)7r//3 in the wide band limit. 

We review the earlier results for the four vertex 
ro,i,2,3('^«oj'^nij'^n2, Wnj), which analytically continued 
to real frequencies and then evaluated at zero frequency. 
Diagrams for this vertex up to third order are shown in 
Fig. 8. In third order of U, the interaction vertex correc- 
tions are given by 



412,3(0,0,0,0) = -3^ ^ GliLO,)G„{uj2)Go{uJ2) 



f/3 



D 



Aa(£)& 
TT J_jj iiUn + l-J. - e 



iAaiiUn + ^)sgn(w„) 
(50) 



-1-3— hni ^ Ga{(^l)Ga{uj2)Go{uj2)Ga{(^l + UJ2 - uj) 



32 



-h3^ lim V Ga{i^i)Ga{uJ2fGo{iUi+uj2-uj), (54) 



where Aa{e) = nV^poie). In the flat wide band hmit, 
/9o(e) is independent of e and D 00, Aa{e) becomes a 
constant. 



which correspond to the Feymann diagrams (b), (c), (d) 
in Fig. 8, respectively, where the Green's function prop- 
agators are all for f/ = 0, as given in Eq. (p3|). Diagrams 
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(b) and (c) can be severed into separate diagrams by 
cutting a pair of lines, one of which is the propagator 
of the a = Majorana fermions. Such diagrams give 
singular contributions to the interaction vertex (parquet 
diagrams). They can be evaluated as follows 



u 

ttA 



^2 2ttT' ^2' 



for T < A. (55) 



While for the second term, wc have to complete the sum- 
mations over the internal frequencies, and then let the 
external frequency go to zero. After that, we obtain the 
following triple integral 



3jj3 roo 

Ic = -7- / deide2dezp{ei)p{e2)p{e-z) 

^ J —00 



2T 



£2(61 + £2 - £3) cosh(|i^) cos}i{^) cosh(ll^) 



(56) 



By changing the variables, we can reduce the integral, 
and then extract the leading singular contribution, which 



and Xj is given by 

Xj -«(.9ji"B) 



2 d{did2) 



(59) 



The first singular contribution to Xj comes from the sec- 
ond order in U, described by the diagram (a) in Fig. 9. 
It is straightforward to evaluate it, and obtain 



(2) 



(60) 



When we consider the fourth order in U, higher loga- 
rithmic terms appear, and they are described by the di- 
agrams (b) and (c) in Fig. 9. According to the previ- 
ous argument, these two diagrams are also parquet ones, 
which will produce the singular contributions to Xj ■ The 
evaluation of the diagram (b) is easier, and we can get 



11, = 



^Eg'(^i) ^EG"(^2)Go(a;2) 
ttA ttA 1 



(61) 



While for the diagram (c), it is more complicated to cal- 
culate its contribution. It is denoted as the following 
term 



dei 



tanh(£i/2T) 



ei 



(57) 



where a high energy cutoff factor A has been introduced 
when calculating the integral. 

The above results have clearly shown that the logarith- 
mic contributions to the interaction vertex corrections in 
the third order perturbation theory cancel exactly. We 
are left diagram (c) in Fig. 8 which contributes a regular 
term only. Therefore, up to third order in U, there are 
no singular vertex corrections. In a renormalized per- 
turbation theory this vertex is multiplied by a wavefunc- 
tion rcnormalization factor which can be shown does have 
ln(T) terms. However, the sign of this term is such that 
in the rcnormalization group equations of the form given 
in jj] the renormalized interaction decreases rather than 
increases so the fixed is stable in the presence of U and 
not unstable as concluded earlier. 

In a similar way we can show that the (ln(T))^ contri- 
butions to the total spin susceptibility Xj that occur to 
order also cancel. This susceptibility is the response 
of the impurity to a field H12 (= gjpBH) coupled to the 
jz component of the impurity states. This coupling cor- 
responds to an extra term in the action of the form, 



iH 



12 



dl {T)d2{T)dT, 



(58) 



/; = C/'(^)'(^)' E G4c.i)Go(a;i)G4^2) 
X Go(a;2)Ga(w3)Go(cJ2 -t^i -t^a)- (62) 



Here we also have to finish the summmations over fre- 
quencies first, and then we find that 



8 VA^ 



(i£ld£2(i£3(i£4p(£l)p(£2)p(e3)/3(£4) 



1 



e3e4(ei + £2 - £3 + £4) 



sinh( ^^+%7-^+'^- ) 



cosh(|y) coshdy) cosh(|y) cosh(|y) 



(63) 



After that, by changing variables we can reduce the in- 
tegral further and identify the leading singular contribu- 
tion, which is 



" VA^ ttA 



^ tanh(£i/2T) 
d€\ - 



£1 



de- 



tanh(£2/2T) 



£2 



ttA ttA 1 



(64) 



Therefore, the squared logarithmic contribution terms in 
the fourth order perturbation for Xj cancel exactly, but 
there are still the logarithmic terms in this order. 



11 



B. Ward Identity 

To derive the Ward identity we exploit the symmetry 
noted earlier that the Hamiltonian is invariant under an 
0(3) orthogonal transformation in the space spanned by 
the 1,2,3 Majorana fermions. We first of all look at the 
effect of using such a transformation to change the vari- 
ables in the action. We look at a particular one parameter 
transformation, corresponding to a rotation through an 
angle 9 in the 1,2 plane about axis 3, 

diir) = cos6I(t)4(t) - sin6'(T)4(T), 

da (r) = sm9{T)d[ (r) + cos6'(T)d2 (r) (65) 

For 6 independent of r, and in the absence of external 
source fields, the partition function is invariant under 
this change of variables. However, when 9 is taken to 
be a local r dependent transformation a new term in 
the action is generated, and there are further new terms 
arising from the couplings to the external source fields, 
as these are also not invariant under the transformation. 
Apart from these extra terms the expression for the gen- 
erating function has the same form in terms of the new 
variables (and the same measure) so these extra contri- 
butions must cancel. The equation for the cancellation 
of these two terms leads to the Ward identity. We work 
to first order only in 9. The extra term to the partition 
function from the coupling to the external source fields 
can be written in the form, 

9iT)irJ^ir){d,iT))~rJ,ir){d^ir))) (66) 

where we have divided by Z in order to express the result 
in the form of an expectation value (with respect to the 
action for 9 — 0). In terms of the Fourier coefficients this 
becomes 

■^'^9{il,n)(r]i{~Uln){di{LJn - ^,n)) - 

m{-^n){d2{uJn-^-rn))^ (67) 

9(Vlm) is the Fourier coefiicient in the expansion of 9{t) 
0{r) = -Y,9{n^)e-^^-\ (68) 

with f7„j — 2TO7r//3 and m an integer, as 9{t) is required 
to satisfy the periodic boundary condition 9{(3) = 9{Q). 
The interaction part of the action is invariant under the 
transformation, even when 9 depends on t, so the other 
extra term generated arises purely from the bilinear term. 
This extra contribution can be written in the form, 

^Y.9{n^rn){dl{~LO,,)d2{^n~^r,^)))Fl2{i^n.^m). (69) 

where Fi2{ujn, ^m) is defined by 



F^p{u^, nra) = ri2„")(L^„) - T^fHuJr. " f^m) (70) 

with Tf'p\LUn) = [G'"^^f^{(^n)]~^ ■ Wc now equate the sum 
of these two terms to zero. We can use the fact that 
the resulting equation holds for arbitrary 9{fl„i), so the 
coefficient of each component must vanish. This leads to 
the equation, 

J2{{dli-i0'n)d2iiVn-^m)m2ii0'r.M 
n' 

= E(^2(c^; - nMi-co'j)) - ini-oj'Mioj'r. - a„))) 

n' 

(71) 

which is the basic Ward identity. The terms in this equa- 
tion all vanish in the absence of the external source terms 
so to derive any useful equations from this for our original 
system without sources we must functionally differentiate 
with respect to 77i(— Wn) and r]2{uJn — flm)- On carrying 
out this differentiation and then equating the external 
source terms to zero we obtain the equation, 

-{{di{-UJri)dl{uJn)) - (c?2(-W„ + fljn)d2iuJn - ^rn))) 

= '^{dl{0Jn)d2{-UJn + flm)dl{-Lj'^)d2{uj'„ - ilm))) 
n' 

xFi2{L^'„,nrn) (72) 

which is an expression which relates the one and two 
particle Green's functions. It can be checked in the case 
U = 0, where it gives a trivial identity. For U ^ we 
can express the two particle Green's function in terms of 
the one particle irreducible four vertex, 

Sll(w„) — £22(^,1 — ilm) = 
n' 

X F^2{Jn. n.ra)Gii{co'„)G22K ~ ^rn) (73) 

which is now a non trivial identity relating the self-energy 
of the one particle Green's function to the irreducible four 
vertex. 

We can obtain another similar identity when the cou- 
pling of the impurity to an external field of the form ( [s^ ) 
is taken into account. Such an interaction induces an off 
diagonal term Ei2(ijJ„) in the self-energy. To relate these 
quantities we have to use Dyson's equation which is now 
in matrix form, 

[r(2.o) _ s]G = I. (74) 

The change in the Green's function SG to first order in 
the applied field is given by 

6G = G[r(2'0)JG(")r(2'0) + s^G, (75) 
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and as the matrices are diagonal in the absence of the 
field this can be simplified to 

Gi2 = GnG22(r(2'0)5G("2^r(2^°) + (76) 

The first term involves 60^2 the off diagonal Green's 
function calculated to first order in H12 but with the 
interaction term set to zero. This is easy to calculate, 
and on substituting the result in the above we find 

Gl2(w«) = Gii(w„)G22(Wn)(i-H"l2 + Sl2(t^n)) (77) 

An alternative expression for Gi2(wn) to first order in 
Hi2 can be deduced from the generating function, by 
taking the appropriate functional derivatives to generate 
Gi2(<^n) and then taking a derivative with respect to the 
external field H12. We obtain the equation, 



Gi2(Ti,r2) 



iH 



12 



/3 Jo 



(r,dl(Tl)d2(T2)rfl(r)d2(T))dT 



(78) 



The expectation value on the right hand side of this equa- 
tion is a two particle Green's function and can be rewrit- 
ten in terms of the irreducible four vertex ri2i2, 

Gi2{iuJn) = Gii{iuJn)G22{ii^n){iHi2 + ^-^^ 



Equating this to our previous result and taking the limit 
Hi2 we find 

3Si2(Wti) 1 V^-p / / / r ' f ' \ 

—^77 = « 2^^1212 C^„,t^„,C^„,CJ„)Gii(w„)G22 W„). 

n' 

(80) 

As T ^ the discrete frequencies w„ can be replaced 
by a continuous variable lo and the summations can be 
replaced integrals so Eq. (pOl) becomes 



9S12H 
idHi2 



duo' 

Ti2l2{uJ,UJ,Uj' ,i^')Gll{uj')G22{io')-Z— ■ 

I'K 



(81) 



We can take the same limit of Eq. (|73|), and if we also 
divide by and take the limit SI — *■ 0, and expression 
similar to (|8l|) will be obtained for the derivative of the 
self-energy with respect to frequency (we can use the fact 
that S]ii((jj) = Ti22{^))- We have to be careful in taking 
the limit O — > because in this limit there is a delta 
function contribution to the integrand at w = 0. This 
arises from the derivative of the sgn functions that orig- 
inate from the imaginary part in Eq. (pO|). The term in 
the integrand of ([73|), 



Fi2(c^',f))Gii(c^')G22(c^'-f^), 

can be rewritten in the form, 

[Gii(c^')-G22(^'-f^)] 



(82) 



[S]n(coO-S]22(^^-»)] 



Gn(c^'-f^)G22(t^'), (83) 



The first term in (g3|) gives the delta function term when 
we take the limit f2 ^ 0. In this limit (RSf) becomes 



Gn(c^')G22(w') + ^(pii(O) + P22(0))<5(co'), (84) 

where PaoSS*) is the spectral density of the Gaa Green's 
function at the Fermi level. Here we have used the fact 
that the Green's functions Gn and G22 are equal and 
the corresponding self-energies are equal but we write 
the result in a symmetrical way. We have also taken the 
flat infinitely wide conduction band limit for simplicity. 
The final expression for the derivative of the self-energy 
is 



9Eii(w) , 95]22(c^) 



iduj 



iduj 



7r(pii(0)+p22(0))ri212(w,W,0,0) 



duj' 

2 / ri212(c^,C^,C^',C^')Gll(c^')G22(^')^, 



(85) 



and hence we find on using (81) the identity. 



d^iijio) ^ 9I]22(^) ^ 2^^i^ 
iduj iduj idHi2 

-7r(pii(0) +p22(0))ri212(w,W,0,0) 



(86) 



These two equations can be used to derive an exact ex- 
pression for the impurity contribution to the susceptibil- 
ity Xj- We calculate Xj using 



Xj 



(ffjMB) 



, d{did2) ^ 
' dHi2 

GnHG22{^) [ 1 



9Ei2(w)\ duj 
idHi2 J 2^' 



(87) 



We can now use (^6|) and substitute the expression for 
the derivative of the off diagonal self-energy. 



Xj 



Gfi(^) 1 



9S]ii(aj)\ duj 



iduj 



Pii(O) / G?i(cc;)ri2i2(w,w,0,0) 



2tt 
duj 
2^ 



where we have used the fact the Green's functions and 
self-energies for the vector channels a = 1,2 are equal. 
The first integral is easy to evaluate because it is the 
derivative of Gii{uj) (w 7^ 0), while for the second in- 
tegral we can use ( |85| ) for = 0, plus the fact that 
ri2i2(0, 0, 0, 0) = due to antisymmetry. We then obtain 
the result. 
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Xj = (5jMb)^/3ii(0) <^ 1 



iduj 



(89) 



which is an exact result for this susceptibility in the limit 
0. 

We can obtain a similar exact result for the specific 
heat coefficient 7 by applying the same line of reasoning 
as is used for the SAM but in terms of the Majorana 
Green's functions. The result for the impurity contribu- 
tion to 7 is 



7 = 



Paa(0)n 



iduj 



(90) 



j=0 



where the sum runs over a = 0, 1, 2, 3 for 14 7^ V/2. For 
Va = V/2 the sum runs only over the vector components 
a = 1,2,3 as the scalar contribution to the free energy 
gives a term in T rather than T^, which contributes to 
the residual entropy term rather than the specific heat 
coefficient. 

From these results we can derive a Wilson ratio (i?) of 
Xjh which is universal in the case K = and Va = V/2. 
The first case corresponds to the 0(4) symmetric model 
in which each Majorana fermion term makes an equal 
contribution so that 



R = 



4^2 fc2 



Xj 



'^{ajmY 1 



(91) 



This result is independent of U in agreement with the 
NRG results. We also know that in the large U limit 
Xj ~^ Xs where Xs is the spin susceptibility, as the isospin 
fiuctuations associated with the impurity are suppressed, 
so this result agrees with the Wilson ratio of the spin 
susceptibility to the specific heat coefficient Xs/l iu the 
large U limit. For 14, = we can derive similar exact 
expression for the spin and charge susceptibilities from 
Ward identities ||l^, because both these quantities are 
conserved when the anomalous hybridization term van- 
ishes. We cannot do so for Va 7^ as these quantities are 
no longer conserved. It is possible to derive some general 
Ward identities for the derivatives of the self-energies of 
all four Majorana fermions when Ki 7^ because the non- 
conserving terms do not involve the two-body interaction 
term jT?! . However, expressions for the spin and charge 
susceptibilities cannot be deduced in terms to the self- 
energies evaluated at zero frequency, as they do not in- 
volve conserved quantities, so that the information given 
by the Ward identities cannot be exploited in the same 
way. 

For Va = V/2 we again find a universal Wilson ratio. 



R = 



47r2fci 



Xj 



3(5j7^b)2 7 



(92) 



The value has changed from the SAM value simply be- 
cause one of the Majorana fermions does not contribute 



to the specific heat coefficient and so the results is | 
rather than |. The corresponding Wilson ratio for the 
spin susceptibility should have the same value in the large 
U hmit. 

We know in this non- Fermi liquid case Vg, ~ V/2 that 
there are singular contributions to both Xj ^^^d 7 as 
T — + 0. This singular temperature dependence is as- 
sociated with the derivatives of the self-energies of the 
vector Majorana fermions evaluated in the limit a; ^ 0, 
and in evaluating the above expressions in this case we 
have to keep a small but finite value for the temperature 
T. We can then define a temperature dependent factor 
Za{T) = z{T) for the vector components a — 1, 2, 3 via 



Zo(T) = 1 - 



iduj 



(93) 



Then the contributions to the specific heat and suscepti- 
bility can regarded as due to free Majorana quasipar- 
ticles each having a renormalized resonance width of 
A(T) = 2;(T)A. In the Wilson ratio the renormaliza- 
tion factor ^(r) cancel so the result is the same as that 
for the non-interacting Majorana fermions. This is simi- 
lar to what happened for the SAM (T4 = 0), except that 
z{T) in that case was independent of T and so the free 
quasiparticle description is the more conventional Fermi 
liquid one. In the calculation of Xj for the SAM the field 
H12 acts only on the up electrons so that the renormal- 
ized interaction U between the quasiparticles (related to 
the vertex roi23(0, 0, 0, 0)) plays no role 

Though we have shown it is possible to interpret the re- 
sults for the non-Fermi liquid low temperature behaviour 
in terms of temperature dependent quasiparticles, we do 
not gain any insight into the nature of the singular terms 
contributing to z(T) other than from the weak coupling 
perturbation theory. Also as l/z{T) 0, which follows 
from the ln(T) of Xjj the quasiparticle weight vanishes 
at T = so in that sense the quasiparticles disappear 
as r — > 0. In the next section we can give a more sat- 
isfactory interpretation of the results by separating out 
the temperature independent contributions to z(T) from 
the singular temperature dependent ones. We can then 
define more conventional quasiparticle excitations. The 
singular scattering of these quasiparticles is not to be in- 
tepretated as a breakdown of the quasiparticle concept. 
It is simply due to the fact that there are real quasipar- 
ticle scattering processes which lead to real singularities 
in the physical properties of this model as T — > 0. 



VII. RENORMALIZED PERTURBATION 
THEORY 

We know that in the case of a Fermi liquid fixed point 
there are quasiparticle excitations from the ground state 
of the interacting system which are in one-to-one cor- 
respondence with the single particle excitations of the 
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non-interacting system. One of us (ACH) showed in an 
earlier paper [T^ that it is possible to make a renormal- 
ized perturbation expansion, valid for all values of U, in 
terms of these quasiparticles about this fixed point for the 
0(4) symmetric Anderson model (with a generalization 
to translationally invariant systems in Ref. ||l^), where 
the exact results for the charge and spin susceptibilities 
at T = are obtained from the lowest order (tadpole) 
diagram and the exact impurity contribution to the con- 
ductivity as T — > from the second order diagrams. For 
[/ = in the 0(3) model we have independent parti- 
cle excitations from the ground state. We now look into 
the possibility of defining quasiparticles and of deriving a 
renormalized perturbation expansion for the 0(3) model 
in the non-Fermi liquid case T4 = V/2, such that the low- 
est order diagrams give exact results in the limit T — > 0. 
We first of all rewrite the self-energy of the Majorana 
fermion retarded Green's functions in the form, 



s„„(c^) = i]('2(^) + si''2(^), 



(94) 



where I]a2(w) is the part of the self-energy at T = 
which gives a contribution to the derivative of I]qq(w) 

which diverges as w — > 0, and si'2(w) is the remaining 
part. Note that as we are using retarded Green's func- 
tions here these are in terms of real frequency variables 
Lu rather than imaginary ones used in the previous sec- 
tion, where we used the thermal Green's functions (they 
are related by G^°^{iujn) — Gaa{i^n)) This latter term we 
write in the form, 

Si'2(c^) = SW(0) + c.S'W(0) + S('r^(c.). (95) 

Substituting this result in the corresponding Green's 
function, and using the fact that from particle-hole sym- 
metry T,al{0) = 0, we find the interacting Green's func- 
tion can be written in the form, 



iAa - tlalii^) - ^aiiio) 



where 



EW(..)=z„Sir^(..), 



(96) 



(97) 



and z„ = 1/(1 - S'ia^(O)). We can rescale the Majo- 
rana field to absorb the factor, and the correspond- 
ing Green's function without this factor will describe the 
corresponding quasiparticle Green's function. Following 
a prescription similar to that used for the 0(4) model we 
can rewrite the Hamiltonian in the form. 



(98) 



where H is the original or "bare" Hamiltonian, as given 
in Eq. (IG), if is a Hamiltonian of the same form where 



all parameters and fields have been replaced by the cor- 
responding quasiparticle ones (indicated by a bar), and 
He consists of the remaining terms, known as counter 
terms. Expressions can be derived for the counter terms, 
which correspond to hybridization with the impurity and 
an on-site interaction term. However, following a mod- 
ified form of the standard renormalization prescription, 
they can be determined by the conditions, 



A 
B 
C 



sL'^2(o) = o, 
s:(?(o) = o, 

C7= (Z0^1Z2Z3)1/%123(0,0,0,0), 



(99) 



which are applied at T = 0. These conditions express the 
fact the quasiparticles and the interaction term is fully 
renormalized with respect to the non-singular contribu- 
tions. One can now develop a renormalized perturba- 
tion expansion for the self-energies si'2(w) and Sa2('^) 
in which one expands both in the interaction term in H 
and all the counter terms in He- These are organized or- 
der by order in powers of U . The condition B expresses 
the requirement that the quasiparticle Green's function 
has weight 1, so the renormalized fields have to be scaled 
appropriately by a change 5za, which is also a counter 
term but it does not occur explicitly in (p^). If, how- 
ever, we use the functional integral formalism then the 
equation corresponding to (^ is 



S = S + Sc 



(100) 



where S is the action, S the renormalized action and 
the counter terms. The counter term contribution in this 
case contains an explicit term in (Jz^, arising from the 
term in (|49|) involving the derivative of the fields with 
respect to r. In this respect the functional integral for- 
mulation of the renormalized perturbation expansion is 
the more natural one and corresponds to the standard 
renormalization procedure for the cj)^ field theory [po| . 

The renormalized perturbation prescription corre- 
sponds simply to a reorganization of the original per- 
turbation theory, no terms have been neglected. Up to 
this point we have been quite general and, if we have 
no singular contributions, the renormalized expansion is 
the same as that for the 0(4) model with 14 = 0, which 
describes an expansion about a Fermi liquid fixed point. 
The prescription here, however, generalizes the expan- 
sion about the Fermi liquid fixed point to the case of any 
14 such that 14 7^ V/2. We now consider the non- Fermi 
liquid case, 14 = V/2. In this case we simplify the no- 
tation and drop the subscripts for a — 1,2,3 and retain 
them only for a = 0. The results for the specific heat 
coefficient 7 and the susceptibility Xj at T = are due 
to the non-interacting quasiparticle system described by 



7rfc| 
2A 



Xj 



ttA 



(101) 
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These give a Wilson xll ratio of |. There are no regu- 
lar corrections to these results at T = (in zero field 
= 0) arising from the renormalized expansion in 
powers of U . This is a consequence of the Ward iden- 
tity we derived in the previous section. We can write the 
wavefunction renormalization factor z(T) used there as 
z{T) = z-Fz('')(r) where z = 1/(1-E'W(0)) and z(^)(T) 
is the remaining part which contains InT terms. The 
result from the non-interacting quasiparticles as defined 
above takes account of all but the singular terms. The 
singular contributions to the susceptibility as we saw in 
the previous section, arise from the diagrams which can 
be severed into separate diagrams by cutting a pair of 
lines, one of which must correspond to the a = prop- 
agator. The results for the second order renormalized 
diagrams for these singular contributions are 



2A \ttAJ \ a J 

{jhm? (US ' 

ttA V ttA 



In 



(102) 
(103) 



These are the only singular terms as T — > 0, as higher 
order weak coupling diagrams to order ln(T), only con- 
tribute to the non-singular renormalized four vertex and 
so are already included in U . There are no counter term 
diagrams to take into account for the singular terms. The 
higher order terms in ln(T) cancel from the same argu- 
ments given in the discussion of the weak coupling theory. 
In the limit C/ — > 0, C7 ^ C/, A ^ A, and we recover the 
weak coupling result. However, the renormalized theory 
is valid for all values of U and so if we take U 3> ttA, in 
which parameter range we can map the model into the a- 
T model then the leading order corrections to the specific 
heat and susceptibility as T ^ are given by the same 
diagrams but now there is only one low energy scale, the 
Kondo temperature Tk, so U ~ ttA. If we put together 
the regular and singular terms we again get a Wilson x/l 
ratio of | independent of U (and U) , as we would expect 
from the Ward identity argument given in the previous 
section. If we compare these results with those derived 
via the Bethe ansatz for the two channel Kondo model ||] 
then the results are the same if U/irA = 1 with A = Tk 
(which also holds for the symmetric Anderson model in 
the Kondo limit 

It is only in the large U limit {U ^ ttA), when the 
isopin fluctuations of the impurity are suppressed, that 
we can identify Xj with the spin susceptibility of the im- 
purity. For U > ttA, we cannot calculate the spin sus- 
ceptibility of the impurity from just the low order terms 
in the renormalized perturbation expansion. This is be- 
cause the spin is not conserved and an expression for 
the spin susceptibility cannot be derived simply from a 
knowledge of the low lying energy levels. The same holds 
true in the NRG approach. The calculation of Xj via the 
NRG is relatively straightforward because the total spin 



plus isospin is conserved and the susceptibility can be 
derived in terms of the energy eigenvalues. The spin sus- 
ceptibility Xs is much more difficult to calculate because 
it requires a knowledge of matrix elements which have to 
be calculated and updated at each step in the NRG calcu- 
lation. There is another way of looking at this difficulty. 
We can calculate the spin susceptibility if we know the 
low lying levels of the system in the presence of a mag- 
netic field which couples only to the impurity spin. In 
principle these levels could be derived from an effective 
Hamiltonian which describes the system near the low en- 
ergy fixed point in the presence of a magnetic field. How- 
ever, there would be another renormalized parameter in 
this Hamiltonian associated with the coupling of the spin 
to the field, which we can interpret as a renormalized g- 
factor. We could calculate the spin susceptibility but the 
g-factor would be unknown. Similar arguments apply to 
the calculation of the charge susceptibility Xc of the im- 
purity. 

There is no problem in calculating Xs and Xc for the 
0(4) model for 14 = as both charge and spin are con- 
served. Exact results are obtained for these quantities 
(at T = 0) from the renormalized perturbation theory 
to first order in JJ |l^. The higher order contributions 
can be shown to cancel as a result of a Ward identity. 
Coleman and Schofield have used the renormalized 
perturbation theory to first order in U to calculate Xs 
and Xc for the model with Va but Va ^ V/2 and 
deduced a value for the Wilson ratio. However, as the 
spin and charge are no longer conserved performing the 
calculation only to first order in U in this case cannot 
be exact. Their result for the Wilson ratio Xs/"f was 
found to be in satisfactory agreement with that given by 
the Bethe ansatz for the anisotropic two channel model 
and to provide an interpolation over the whole parameter 
range (though it differs by a factor of two from the Bethe 
ansatz result in the limit (Ji — J2) ^ Ji)- 

We can give an alternative argument to that of Cole- 
man and Schofield using renormalized perturbation the- 
ory to calculate both x and 7 in the limit Aq ^ A, 
U ^ ttA. In this limit logarithmic terms in Aq ^ A 
appear associated with the same diagrams that give the 
logarithmic in T contributions for Aq = 0. We can reor- 
ganise the renormalized perturbation theory for Aq <^ A 
so that Y^al (i^) refers to the diagrams that give logarith- 
mic terms in Aq as Aq 0. We calculate asymptotically 
the leading order contributions to the specific heat and 
susceptibility in a similar way to the above calculation. 
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U 
ttA 



ln(Ao/A) 
ttA 



ln(Ao/A) 
ttA 



Ao 



(104) 



(105) 
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In the specific heat term the most singular contribu- 
tion as Aq ^ is the one proportional to Aq"^ and so 
dominates in the Wilson ratio in this limit. We found 
earlier that the results in the Kondo limit U ^ ttA 
iXj ~^ Xs) correspond to the two channel Kondo model 
for U/ttA = 1 so we can use this value in the coefficient 
of the leading singular term. The Wilson ratio then gives 
-8(Ao/A)ln(Ao/A) as (Aq/A) ^ in complete agree- 
ment with the Bethe ansatz result for the anisotropic two 
channel model in this limit |^ where we identify ( Aq / A) 
as the ratio of the renormalized energy scales Ta/Tj. The 



Wilson ratio based on (104) and (105) will also be cor- 



rect for the 0(4) model which corresponds to the limit 
Aq A. The logarithmic terms then vanish and we 
recover a value of 2 obtained earlier {U /ttA = 1 in this 
limit also but this factor plays no role in the result as 
it multiplies the log terms only). Hence this result is in 
complete agreement with those derived from the Bethe 
ansatz in both limits, and also provides an interpolation 
between them. 



VIII. CONCLUSIONS 

The two models we have been considering here were 
both introduced as a way of gaining insight into the na- 
ture of the low energy fixed point for the isotropic and 
anisotropic two channel Kondo models. Our results con- 
firm that there are very strong similarities in the results 
for these compactified models (in the localized or Kondo 
limit) and the two channel Kondo model, and that the 
low temperature behaviour of the compactified models 
can be given a simple interpretation in terms an effective 
Majorana fermion model at the low energy fixed point. 
A consistent picture of the isotropic and anisotropic fixed 
points emerges from the perturbational, numerical renor- 
malization group, conformal field theory, and renormal- 
ized perturbation theory approaches. There are, how- 
ever, some important differences which have not been 
brought out in previous work on these models. 

In the NRG results for the isotropic two channel model 
pl| , ^ it has been found that as the fixed point is ap- 
proached that there is no change of the spectrum as 
the number sites N included in each of the two chains 
changes from iV to + 1. This is in contrast to the 
equivalent results for the single channel Kondo model 
where the spectrum of low lying many-body levels varies 
as to whether N is odd or even, and where there are 
two equivalent fixed points corresponding to a transfor- 
mation N N + 2, one for N odd and the other for 
N even. Superficially it appears that there is also no 
change in the spectrum in the NRG results for the com- 
pactified models for the isotropic case as iV ^ -I- 1 
as the fixed point is approached. However this is not 
the case. The degeneracies of the excitations depend on 
whether A^ is odd or even. The levels and the degen- 



eracies can both be explained on the basis of free vector 
and scalar Majorana fermion chains with different bound- 
ary conditions as described in Sec. IV. When A is even 
the many-body excitations correspond to those given in 
Table VII, which are built up from those of a scalar Ma- 
jorana fermion chain with periodic boundary conditions 
and a vector one with antiperiodic boundary conditions. 
When A^ is odd they correspond to those given in Table 
VIII where the boundary conditions on the scalar and 
vector components have been interchanged. As in the 
one channel Kondo case, these are alternative situations 
and cannot be superimposed on each other. However, to 
reproduce the excitation spectrum as found in the NRG 
calculations for the two channel model these alternative 
situations (refered to as sector I and sector II) have to 
be combined. The degeneracies of the many-body ex- 
citations of the two channel model can be explained by 
including the extra uncoupled degrees of freedom in the 
form of an extra non-interacting fermion chain. As this 
is uncoupled from the impurity like the scalar Majorana 
fermion chain it is taken to have the same boundary con- 
ditions as the scalar part. With this model both the 
many-body energies and the degeneracies of the NRG 
spectrum at the fixed point can be reproduced. Only the 
ground state degeneracy, which only acts as a multiplica- 
tive factor for the degeneracies of the excited states, has 
to be reduced as compared to the compactified model. 

The single particle levels of the combined chains with 
the appropriate boundary conditions are shown in Fig. 
10 (for comparison with the spectra of the compactified 
model see Fig. 11). The resulting many-body energies 
with the degeneracies are given in Table IX for sector I 
and Table X for sector II. They agree with those given in 
the NRG calculations of Cragg et al. and Pang and 
Cox an d the conformal field theory results of Affleck 
et al. |9|Jl5|] . This spectrum has also been obtained for the 
two channel Kondo model from single particle excitations 
corresponding to different boundary conditions by Ye |2^ ] 
using the bosonization approach. Ye in a recent paper has 
also constructed the excitation spectrum for the fixed 
point of the compactified models in a similar way [^ . 
His results differ from ours as he has the two sectors, 
which we identify with A^ odd and A^ even, combined as 
in the two channel case. This is not consistent with the 
NRG results. 

Though Affleck and Ludwig, in their conformal field 
theory for the two channel Kondo model, could explain 
the excitation spectrum found in the NRG calculations 
at the fixed point, they could not derive the many-body 
states from a single particle picture with modified bound- 
ary conditions. Nor could they obtain an explicit form for 
the leading irrelevant operator at the fixed point in terms 
of local operators. We have shown, however, that both 
these can be achieved within the conformal field theory 
approach for the cr-r model. In Sec. Ill we found the 
same two separate sectors of excitations, built up from 
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free Majorana chains with different boundary conditions 
as in the NRG approach, and the same local operator as 
in the NRG was found for the leading irrelevant interac- 
tions using very similar arguments to those used for the 
single channel Kondo model. With these insights from 
the CFT of the cr-r model it might prove possible to get a 
deeper understanding of the two channel conformal field 
theory, and to understand why the sectors correspond- 
ing to Tables VII and VIII have to be combined for that 
system. 

At this point it would seem appropriate to compare 
and contrast our results with those of previous papers on 
these models. We agree with the work of Coleman et al. 

that the fixed point of the isotropic cr-r model does 
correspond to strong coupling and we agree with the form 
of the leading irrelevant interaction. However we find 
that D/J is not a good expansion parameter as {D/JY 
corresponds to U /ttK which is of order unity. We have 
also been able to resolve the problems of the higher order 
logarithmic terms that were found in our earlier work 
on this model (G-M.Z., A.C.H.) ||] in the weak couphng 
expansion for the 0(3) symmetric Anderson model (Vq = 
0). Explicit cancellation of these higher order log terms 
to the vertex roi23(0, 0, 0, 0) and the susceptibility Xj has 
been demonstrated up to fourth order, and we believe 
this can be generalized to all orders in U . This means 
that the marginal Fermi liquid fixed point is marginally 
stable rather than marginally unstable. 

As the nature of the fixed point found in the NRG 
calculations is always independent of [/, we have found 
it natural to describe the fixed point as a renormalized 
form of the Anderson model. This is similar to the case 
of the standard Anderson impurity model where the fixed 
point always corresponds to a Fermi liquid fixed point, 
whatever the value of J7, and the low temperature be- 
haviour can be described by a renormalized version of 
the same model. In the strong correlation regime the pa- 
rameters of the renormalized model are not independent. 
They depend on the single energy scale set by the Kondo 
temperature Tk. We find similar results for the marginal 
Fermi liquid fixed point of the 0(3) model. The main dif- 
ference with the standard Anderson model is that there 
is a singular scattering mechanism of the quasiparticles 
which leads to real singularities in some of the proper- 
ties of the model, such as the susceptibility, as T ^ 0. 
This can be described by a modified form of the renor- 
malized perturbation theory that was developed for the 
0(4) model [|l8|, in which the low temperature properties 
are given exactly from the diagrammatic expansion up to 
second order. 

For the 0(3) Anderson model with < Vo < ^ hi the 
strong correlation regime the renormalized model has two 
energy scales rather than one, as does the corresponding 
cr-r model and the anisotropic two channel model 
Apart from the fact that there are two independent en- 
ergy scales the fixed point behaviour is similar to that 



for the 0(4) model and could be described as a Fermi 
liquid fixed point (as in |25j). As the spin and charge are 
not conserved for this model we found that the first order 
renormalized perturbation theory for the spin and charge 
susceptibilities, as used by Coleman and Schofield ||], is 
not exact as it is in the case of the 0(4) model (where 
they are conserved). However, we found that it is pos- 
sible within this approach to calculate the Wilson ratio 
asymptotically in the limit Vg — + which gives a result 
in complete agreement with that from the Bethe ansatz 
for the anisotropic Kondo model in the same limit. 

The 0(3) symmetric Anderson model is of interest 
apart from its relation to the two channel Kondo model as 
it displays a marginal Fermi liquid fixed point for Vq = 0. 
The resistivity for this model does not correspond to the 
resistivity or the two channel Kondo model and conse- 
quently has a different temperature dependence. In the 
lowest order perturbation theory it was shown earlier Q| 
that to lowest order in U it is linear in T term, as in 
the marginal Fermi liquid theory. The dynamics of this 
model over the full parameter range is presently being 
calculated using the NRG . 
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TABLE L Many-body excitation energies and their corre- 
sponding degeneracies for a single Majorana fermion chain 
with periodic boundary conditions. The single particle spec- 
trum is given in Fig. 1(a). 
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TABLE II. Many-body excitation energies and their cor- 
responding degeneracies for a single Majorana fermion chain 
with antiperiodic boundary conditions. The single particle 
spectrum is given in Fig. 1(b). 



Ee^/{nVF/l) 




dg 


total dg 








1 


1 


1/2 


£l/2 


1 


1 


3/2 


£3/2 


1 


1 


2 


£1/2 + £3/2 


1 


1 


5/2 


£5/2 


1 


1 


3 


£1/2 + £5/2 


1 


1 


7/2 


£7/2 


1 


1 


4 


£1/2 + £7/2 
£3/2 + £5/2 


1 
1 


2 


TABLE HI. Many-body excitation energies and their cor- 
responding degeneracies for a vector Majorana fermion chain 
with periodic boundary conditions. The single particle spec- 
trum is given in Fig. 2(a). 
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TABLE IV. Many-body excitation energies and their cor- 
responding degeneracies for a vector Majorana fermion chain 
with antiperiodic boundary conditions. The single particle 
spectrum is given in Fig. 2(b). 


E^^/{nVF/l) 




dg 


total dg 








1 


1 


1/2 


£1/2 


3 


3 


1 


2£l/2 


3 


3 



3/2 £3/2 3 

3£i^2 1 4_ 

_2 £1/2 + £3/2 9 9_ 

5/2 £5/2 3 

2£i/2 + £3/2 9 12 

3 £1/2 + £5/2 9 
3£i/2 + £3/2 3 12 

7/2 £7/2 3 

2£i/2 + £5/2 9 21 

£1/2 + 2£3/2 9 



19 



TABLE V. Finite-size spectrum of the non-interacting case 
for PBC (with degenerate ground state) . Compared with the 
APBC case, the ground state energy shifts by 1/4 unit, dg is 
the degeneracy of the energy level. 
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TABLE VI. 


Finite-size spectrum of the non-interacting 


case for APBC (with 


non-degenerate ground state). 
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TABLE VII. Finite-size spectrum for J = J* correspond- 
ing to the vector part with APBC and scalar part with PBC. 
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TABLE VIII. Finite-size spectrum for J = J* correspond- 
ing to the vector part with PBC and scalar part with APBC. 
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TABLE IX. Many-body excitation energies and their cor- 
responding degeneracies for sector I of the two channel Kondo 
model corresponding to the single particle spectrum of Fig. 
10(a). 
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TABLE X. 


Many-body excitation energies 


and their corre- 


sponding deg 


;eneracies for sector II of the two channel Kondo 


model corresponding to the single 


particle spectrum of Fig. 


10(b). The energies are measured relative to the ground state 


of sector I. 
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FIG. 4. The picture of the strong-coupling fixed point 
Hamiltonian in the lattice version. The impurity spin com- 
bines one electron at the nearest site to form a singlet bound 
state, and both the impurity spin and the nearest site of the 
vector part thus drop out of the problem. 



FIG. 1. Single particle spectrum of a single component 
Majorana fermion chain; (a) for periodic boundary conditions, 
(b) for antiperiodic boundary conditions. Energies are give in 
units of -kvf/1- 
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FIG. 5. Flow diagram for U = 0.001428, V = 0.01414 and 
Va = V/2. Solid and dashed lines belong to jtot = 0.5 and 
itot = 1.5, respectively. The system flows to the non- Fermi 
liquid fixed point. 



FIG. 2. Single particle spectrum of a vector Majorana 
fermion chain; (a) for periodic boundary conditions, (b) for 
antiperiodic boundary conditions. Energies are give in units 

of TTVf/1- 
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FIG. 3. The separation of the model in the lattice form. 
The upper chain represents the vector part of the conduc- 
tion electrons consisting of ipi, V'2, and ips, while the lower 
chain represents the scalar part tpo, which decouples from the 
impurity spin. 
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FIG. 6. A-dependence of the many-body energies -Bex. The 
ground state (i5ox = 0) has four-fold degeneracy. Degeneracies 
of the excited states plotted here are dg — 4 (dashed line) and 
dg — 12 (solid line). The points where the i?ox(A) cross the 
vertical line at A = 2.5 correspond to the energies obtained 
from the NRG calculation. 
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FIG. 7. Structure of the non-Fermi liquid fixed point. U 
couples the vector part at site -1 to an additional Majo- 
rana-fermion which itself is not coupled to the scalar chain. 
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FIG. 9. The total spin susceptibility in the second and 
fourth order perturbation theory, (a) describes the second 
order correction, giving rise to a logarithmic term; (b) and 
(c) are the leading order singular contributions in the fourth 
order perturbation theory, producing the squared logarithmic 
corrections. 
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FIG. 10. Single particle spectra for the two sectors present 
at the non-Fermi liquid fixed point of the two channel Kondo 
model, (a) corresponds to sector I and (b) to sector II. The 
energies are given in units of 




FIG. 8. The interaction vertex function in the first and 
third order perturbation theory, (a) is the original interaction; 

(b) , (c), and (d) are the third order corrections, (b) and 

(c) are parquet diagrams, producing the leading logarithmic 
terms. 
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FIG. 11. Single particle spectra of the compactified models. 

(a) for PBC in the vector part and APBC in the scalar part 
(corresponding to N odd) and (b) for APBC in the vector 
part and PBC in the scalar part (corresponding to N even). 
The energies are given in units of 
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